Fluctuation-Induced Heat Release from Temperature-Quenched Nuclear Spins near a 

Quantum Critical Point 
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At a quantum critical point (QCP) — a zero-temperature singularity in which a line of continuous 
phase transition terminates — quantum fluctuations diverge in space and time, leading to exotic 
phenomena that can be observed at non-zero temperatures. Using a quantum antiferromagnet, we 
present calorimetric evidence that nuclear spins frozen in a high-temperature nonequilibrium state 
by temperature quenching are annealed by quantum fluctuations near the QCP. This phenomenon, 
with readily detectable heat release from the nuclear spins as they are annealed, serves as an excellent 
marker of a quantum critical region around the QCP and provides a probe of the dynamics of the 
divergent quantum fluctuations. 

PACS numbers: 75.30.Kz, 75.40.-s, 75.50.Ee, 76.60.Es 



Quantum criticality, the divergence of spatial and tem- 
poral extents of quantum fluctuations at a continuous 
phase transition at zero temperature, is the driving mech- 
anism behind a variety of novel phenomena that challenge 
conventional theoretical approaches to collective behav- 
ior of many-body systems and may also lead to technolog- 
ical applications 0, Q • This criticality is responsible for 
the breakdown of Fermi-liquid behavior in heavy-fermion 
metals 0, H[ and for the emergence of exotic states of 
matter such as unconventional superconductors @,@], ne- 
matic electron fluids [?J in GaAs/AlGaAs heterojunctions 
H and Sr3Ru207 0], and yet-to-be-identified 'reentrant 
hidden-order' states in URu2Si2 [ltj. 

Unlike in classical criticality, which is the divergence 
of order-parameter fluctuations in space at a continu- 
ous phase transition occurring at a non-zero tempera- 
ture [Til ] . dynamic and static pro per ties are inextrica- 
bly linked in quantum criticality [H, [13]. A spectac- 
ular demonstration of this linkage is the quantum an- 
nealing — the quantum-fluctuation-driven relaxation of 
quenched disorder — of an Ising spin glass near its quan- 
tum critical point (QCP), a phenomenon with broad im- 



m m 



plications in efficient optimization algorithms 
Here we report the finding of another relaxation phe- 
nomenon near a QCP: quantum- fluctuation-driven re- 
lease of heat from nuclear spins. This finding opens up 
a unique avenue to investigate the dynamics of quantum 
fluctuations that underlie quantum criticality. 

For this study, we have chosen the inorganic coordina- 
tion compound Cr(diethylenetriamine)(0 2 )2-H 2 0, here- 
after referred to as Cr(dien), because of its large num- 
ber of hydrogen nuclear spins \v\ . This quasi- two- 
dimensional quantum magnet has a magnetic-field-tuned 
QCP at the field H c = 12.3 T, where a highly polar- 
ized antiferromagnetic phase gives way to a field-induced 
ferromagnetic phase. The Cr(diethylenetriamine)(02)2 
molecule, the key building block of this material, is an 



oblate, elongated disk-shaped molecule illustrated in Fig. 
[IJa). The magnetic ion is Cr(IV), located on the mirror- 
symmetry plane of the molecule. In the monoclinic crys- 
tal structure of Cr(dien), the S = 1 spins of Cr(IV) form 
a square lattice along the crystallographic ac plane, with 
an exchange energy J of 2.71-2.88 K. The spins order an- 
tiferromagnetically at Tn = 2.55 K in zero magnetic field 
[l7j |. Application of a high magnetic field depresses the 
ordering temperature, pushing it to zero at the critical 
field H c . We work near this QCP. 

The experiment is done in a relaxation calorimeter 
[HI, [HI]- A 1.02 mg single crystal of Cr(dien) is placed 
on the sample platform, weakly linked to the thermal 
reservoir through the leads for the thermometer and 
heater. This geometry allows the sample temperature to 
be raised with the heater with ease and to drop rapidly 
to the reservoir temperature when the heater is turned 
off. It also enables the detection of spontaneously re- 
leased heat in the sample as an observable temperature 
difference between the sample and the reservoir. 

The procedure of the experiment is illustrated in Fig. 
QJb). First, the sample is heated from temperature To 
(point Pi or P[ in the figure) of the thermal reservoir to 
temperature T q (point P2 or P%) ranging from 266 mK to 
1.52 K, in magnetic field H q applied perpendicular to the 
ac plane of the crystal. After 1.4 min at T q , the sample 
is rapidly cooled back to T in 1.8 s by turning off the 
heater. This temperature quenching leaves the hydro- 
gen nuclear spins frozen in a nonequilibrium high-energy 
state corresponding to T q . Subsequently, the magnetic 
field is swept at 0.2 T/min or 0.1 T/min through the crit- 
ical field H c . The field sweeps are made at four different 
Tq ranging from 96 mK to 261 mK. 

As the swept field approaches H c , heat is released in 
the sample, manifesting itself as a pronounced peak in 
the temperature difference AT between the sample and 
the thermal reservoir, as shown in Figs. [He) andQJd). 
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FIG. 1: (color online). (a) Pair of H 2 and 
Cr(diethylenetriamine)(02)2 molecules, the basic unit of the 
Cr(dien) crystal. The Cr(IV) ion is the largest sphere to the 
right, bonded to four oxygen and three nitrogen atoms, (b) 
Procedure of the experiment (see text). Broken line, termi- 
nating in a quantum critical point at T = 0, is the boundary 
between the highly polarized paramagnetic phase and anti- 
ferromagnetic phase, (c) Temperature difference AT between 
the sample and the thermal reservoir as a function of the 
magnetic field during field sweeps at 0.2 T/min. The thermal 
reservoir is held at 181 mK. Temperature T q from which the 
sample has been quenched ranges from 266 mK to 1.52 K. (d) 
Evolution of AT when the thermal reservoir is held at 96 mK, 
as the field is swept at 0.1 T/min. T q ranges from 300 mK to 
797 mK. In (c) and (d), the peaks to the right are observed 
during downward field sweeps and those to the left during 
upward sweeps. The fields indicated in the legends are H q . 



The heat release occurs only during the first field sweep 
after the temperature quenching of the sample from T q , 
not during subsequent sweeps [20 ]. 

The amount of released heat Q is obtained from the 
data via 



Q = / nATdH/H, 



(1) 



where k is the thermal conductance of the weak link be- 
tween the sample and the thermal reservoir, and H the 
field-sweep rate. As shown in Figs. Eta) and[2^b), Q de- 
pends on T q , H q , and To. These dependences indicate 
unambiguously that the heat is indeed released from the 
hydrogen nuclear spins, as the following analysis reveals. 

When hydrogen nuclear spins in magnetic field H q are 
initially frozen in a nonequilibrium high-energy state de- 
termined by temperature T q — while the lattice cools 
rapidly to temperature To — and subsequently equili- 
brate with the lattice at another field H , the amount of 




FIG. 2: (color online). Amount of heat released at (a) 181 mK 
and (b) 96 mK as a function of T q , the temperature from 
which the sample has been rapidly quenched. In (a), the 
sample has been quenched at H q = 6T (o), 9T (♦), 11.5 T 
(▼), 13.5 T (a), 16 T (•), and 18 T (■). The lines represent 
Eq.[2]with n H = 10: dash-dotted line for H q = 11.5 T, dotted 
line 13.5 T, broken line 16 T, and solid line 18 T. In (b), the 
sample has been quenched at 10 T (•) and 13.5 T (■). The 
lines representing Eq. [5] are with ng = 10 for H q = 10 T 
(broken line) and with uh = 15, of which five are assumed to 
freeze only at 132 mK, for 13.5 T (solid line.) 



heat released by them is 



Q = jinnR 



1(1 + 1) 



H 

n 



H. 
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Here nn is the number of hydrogen nuclear spins (per 
formula unit) that participate in heat release, n sample's 
mole number, R the gas constant, h the Planck constant, 
ks the Boltzmann constant, and 7 and 1—1/2 the gy- 
romagnetic ratio and spin of the hydrogen nucleus. The 
equation takes into account that the nuclear-spin tem- 
perature just before the heat release is HT q /H q instead 
of T 9 , as a result of the field sweep from H q to H being 
adiabatic for the nuclear spins except very near H. 

For To = 181 mK, good agreement between experiment 
and Eq. [2] is obtained with uh — 10, as shown in Fig. 
[2{a), except for a few points for which H q is 18 T or 16 T 
|2ll]. At T = 96 mK, the peak that appears at 12.76 T 
in Fig. HJd) during downward field sweeps from 13.5 T 
is very sharp while AT > 36 mK, i.e., while the sample 
temperature T + AT is higher than 132 mK. This indi- 
cates that hydrogen nuclear spins whose relaxation times 
r are very short when T > 132 mK now participate in 
heat release. When To = 181 mK, they have evidently 
reached thermal equilibrium with the lattice before each 
field sweep starts and therefore do not participate in heat 
release. As shown in Fig. [2jb), good agreement between 
Eq. [2]and the T = 96 mK data for H q = 13.5 T is ob- 
tained with nn = 15, of which five whose r are very short 
at To > 132 mK are assumed to freeze at a nonequilib- 
rium high-energy state at 132 mK instead of T q . When 
quenched at H q = 10 T, in the antifcrromagnetic phase, 
riff = 10 gives better agreement with the data, suggesting 
that those five hydrogen nuclear spins do not participate 
in heat release even at Tq = 96 mK. 
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FIG. 3: (color online). Magnetic fields and temperatures 
of the peaks in the AT curves (squares), marking a quan- 
tum critical region delimited by two straight lines through 
the data points. Solid squares are for peaks at a sweep rate 
of O.lT/min, open squares at 0.2T/min. Circles represent 
the phase boundary — detected by the magnetocaloric effect 
— between the antiferromagnetic (AF) and highly polarized 
paramagnetic (P) phases, with the broken line from a power- 
law fit of data points up to 0.84 K. At zero temperature, the 
P phase turns into a field-induced ferromagnetic phase. 



Cr(dien) contains fifteen hydrogen atoms per formula 
unit, as shown in Fig. [IJa). Among the thirteen in the 
Cr(diethylenetriamine)(02)2 molecule, the five bonded 
to nitrogens are closer to the Cr(IV) ion than eight that 
are bonded to carbons [22j . It is very likely that the five 
hydrogens with short nuclear-spin r are those bonded to 
the nitrogens and thus experience stronger fluctuating 
dipolar field and transferred hyperfine field of the Cr(IV) 
ion, whereas the ten hydrogens with longer r are those 
bonded to the carbons and the two in the water molecule. 

The positions of the heat-release peaks [23[ are shown 
in Fig. [3] along with the phase boundary between the 
highly polarized antiferromagnetic phase and similarly 
highly polarized paramagnetic phase, which turns into 
the field-induced ferromagnetic phase at zero tempera- 
ture, of Cr(dien) near its QCP at H c . This diagram 
suggests that, in the zero-temperature limit, the loci of 
the peaks converge on the QCP. These loci delimit the 
region in which r is shorter than the timescale of the 
experiment, a quantum critical region. 

How does r vary within and near this region? To an- 
swer this question, we stop the field sweep at field H 
and record the subsequent relaxation of the sample tem- 
perature toward To [24;]. The measurements are carried 
out at To = 181 mK after the sample has been quenched 
from T q = 796 mK at H q = 13.5 T or from T q = 761 mK 
at H q — 10 T. To remove uninteresting contributions 
coming from the magnetocaloric effect and eddy-current 
heating, the measurements are repeated without tem- 
perature quenching, and the control data thus obtained 
are subtracted from the original data. The relaxation 
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FIG. 4: (color online). Relaxation of the temperature dif- 
ference between the sample and the thermal reservoir, as 
the magnetic field is brought close to H c after temperature 
quenching to To = 181 mK. (a) Relaxation curve after a down- 
ward field sweep from H q = 13.5 T toward H c at O.lT/min 
has been stopped at 12.7 T. The sample has been prepared by 
quenching it from T q — 796 mK at H q . The line is a fit to an 
exponential relaxation with r = 42.6 s. (b) Relaxation rate 
1/tvs. the field at which field sweep has been stopped. Solid 
symbols are data from curves showing exponential relaxation, 
and open symbols from curves showing stretched-exponential 
relaxation. Lines are guides to the eye. (c) Amount of heat 
released vs. the distance from H c . Horizontal lines indicate 
the amounts of heat released during complete field sweeps 
through H c . In (b) and (c), circles and broken lines are for 
= 10 T, and squares and solid lines for 
13.5 T. 
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is exponential at some H as shown in Fig. 0Ja) but 
shows stretched-exponential behavior at other H, with 
the stretching exponent ranging from 0.61 to 0.78. The 
relaxation rate 1/r diverges as the field approaches H c , 
as shown in Fig.QJb), indicating that diverging quantum 
fluctuations of the Cr(IV) spins near the QCP drive the 
relaxation of the nuclear spins. The divergence is asym- 
metric around H c , faster for H > H c than for H < H c . 
This asymmetry is also seen in Fig.[3]as a wider quantum 
critical region for H > H c than for H < H c . 

The amount of heat released during the relaxation-time 
measurements is shown as a function of \H — H c \ in Fig. 
IDJc) . The result suggests that the ten hydrogen spins that 
release heat at this temperature contain two groups, each 
comprising four to six spins per formula unit. Above H c , 
one group relaxes at \H — H c \ of about 0.8 T, whereas 
the other group relaxes at fields closer to H c . Similarly, 
below H c , the first group relaxes at \H — H c \ of about 
0.6 T, whereas the second group relaxes at fields closer 
to H c . The result is consistent with the molecular struc- 
ture of Cr(diethylenetriamine)(02)2 shown in Fig. [Ha): 
the eight hydrogens bonded to the four carbons — with 
longer nuclear-spin r — fall into two groups, each consist- 
ing of four hydrogens, with distinct ranges of distances 
from the Cr(IV) 0. 

In conclusion, our results provide unambiguous evi- 
dence that temperature quenching of Cr(dien) leaves the 
hydrogen nuclear spins frozen in a nonequilibrium high- 
energy state and, as the magnetic field is then brought 
close to the QCP, quantum fluctuations of the Cr(IV) 
ionic spins quickly anneal them to reach thermal equi- 
librium with the lattice. These results also imply that 
the quantum-fluctuation-driven heat release from nuclear 
spins is a generic phenomenon to be found near a variety 
of QCP. Because of the inextricable link between dynamic 
and static properties in quantum criticality, quantum- 
critical systems are predicted to exhibit interesting, non- 
trivial relaxation phenomena during and after a sweep of 
a control parameter such as magnetic field and pressure 
through the QCP [26f and also after temperature quench- 
ing near the QCP [2jJ. Our results warn, however, that 
the response of nuclear spins — which are nearly ubiq- 
uitous — to those changing parameters and to quantum 
fluctuations must be carefully taken into account in real 
solids. At the same time, heat release from nuclear spins 
promises to be a useful probe for the dynamics of quan- 
tum fluctuations that underlie quantum criticality in a 
variety of systems. 

We thank S. Nellutla and A. Kumar for contributions 
to an early stage of this work and K. Ingersent, P. Ku- 
mar, Y. Lee, D. L. Maslov, and R. Saha for useful dis- 
cussions. Thanks are also due to J.-H. Park, T. P. Mur- 
phy, and G. E. Jones for assistance and R. J. Clark for 
drawing a figure. This work was supported in part by 
NSF grant NIRT-DMR-0506946 and by the University of 
Florida Physics REU program under NSF grant DMR- 



0552726. The experiment was performed at the National 
High Magnetic Field Laboratory, which is supported by 
NSF Cooperative Agreement DMR-0654118, by the State 
of Florida, and by the Department of Energy. 



[1] J. A. Hertz, Phys. Rev. B 14, 1165 (1976). 

[2] S. Sachdev, Quantum Phase Transitions (Cambridge 

Univ. Press, Cambridge, 1999). 
[3] H. v. Lohneysen, A. Rosch, M. Vojta, and P. Wolfle, Rev. 

Mod. Phys. 79, 1015 (2007). 
[4] P. Gegenwart, Q. Si, and F. Steglich, Nat. Phys. 4, 186 

(2008). 

[5] N. D. Mathur et al, Nature 394, 39 (1998). 
[6] H. Q. Yuan et al, Science 302, 2104 (2003). 
[7] E. Fradkin and S. A. Kivelson, Phys. Rev. B 59, 8065 
(1999). 

[8] M. P. Lilly et al, Phys. Rev. Lett. 82, 394 (1999). 
[9] R. A. Borzi et al, Science 315, 214 (2007). 

[10] Y. S. Oh et al, Phys. Rev. Lett. 98, 016401 (2007). 

[11] C. Domb, The Critical Point: a Historical Introduction 
to the Modern Theory of Critical Phenomena (Taylor & 
Francis, London, 1996). 

[12] S. L. Sondhi, S. M. Girvin, J. P. Carini, and D. Shahar, 
Rev. Mod. Phys. 69, 315 (1997). 

[13] Q. Si, Adv. Solid State Phys. 44, 253 (2004). 

[14] J. Brooke, D. Bitko, T. F. Rosenbaum, and G. Aeppli, 
Science 284, 779 (1999). 

[15] G. E. Santoro, R. Martonak, E. Tosatti, and R. Car, 
Science 295, 2427 (2002). 

[16] D. A. House and C. S. Garner, Nature 208, 776 (1965). 

[17] C. M. Ramsey et al, Chem. Mater. 15, 92 (2003). 

[18] H. Tsujii et al, Physica B 329-333, 1638 (2003). 

[19] The thermometers of the calorimeter have been self- 
calibrated in magnetic fields, by measuring the heat ca- 
pacities of standard materials in situ. In the field range 
of the present experiment, calibration uncertainties are 
about 1.0% at 0.8 K and about 1.8% at 0.1 K. 

[20] AT during the subsequent sweeps, exhibiting only the 
magnetocaloric effect of at most 1 mK and eddy-current 
heating by 8 mK (T = 181 mK) and 3 mK (To = 96 mK), 
have been subtracted from the data shown in Figs. 1 and 
2. 

[21] For those points, 1.4 min of waiting is probably insuffi- 
cient to thermalize the hydrogen nuclear spins. 

[22] The five hydrogens bonded to the nitrogens are located 
at 2.40-2.52 A from the Cr(IV), whereas the eight that 
are bonded to the carbons fall into two groups of four, 
each at 3.08-3.29 A and 3.77-3.79 A from the Cr(IV). The 
water hydrogens are at 3.29 and 3.48 A from the closest 
Cr(IV) ions, which are in Cr(diethylenetriamine)(02)2 
molecules not shown in Fig. []J a). See supporting online 
information for Ref. [Th] . 

[23] The peak position depends slightly on the quenching con- 
ditions T q and H q . For the data shown in Fig. [3] in 
the highly polarized paramagnetic phase, H q = 13.5 T 
and T q = 0.80 K except for the 227 mK point for which 
T q — 1.04 K. For those in the antiferromagnetically or- 
dered phase, H q — 9-10.5 T and T q — 0.76 K except for 
the 226 mK point (T q = LOOK) and the 201 mK point 
taken at 0.2T/min (T q = 0.50 K). 



5 



[24] t measured in this manner is approximately T\ + C„/k 
[25j ] . where T\ is the nuclear spin-lattice relaxation 
time and C n the nuclear heat capacity. At 181 mK, 
Cn/re =11.1 s at 12.1 T. 

[25] B. Andraka and Y. Takano, Rev. Sci. Instrum. 67, 4256 



(1996). 

[26] D. Patane et al, Phys. Rev. Lett. 101, 175701 (2008). 
[27] D. Patane, A. Silva, F. Sols, and L. Amico, Phys. Rev. 
Lett. 102, 245701 (2009). 



